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Abstract 

We consider an integrable model describing an Anderson-like impurity coupled to an open t— 
J chain. Both the hybridization (i.e. its coupling to bulk chain) and the local spectrum can be 
controlled without breaking the integr ability of the model. As the hybridization is varied, holon 
and spinon bound states appear in the many body ground state. Based on the exact solution we 
study the state of the impurity and its contribution to thermodynamic quantities as a function 
of an applied magnetic field. Kondo behaviour in the magnetic response of the impurity can be 
observed provided that its parameters have been adjusted properly to the energy scales of the holon 
and spinon excitations of the one-dimensional bulk. 



I. INTRODUCTION 



The possibility of controlled embedding of quantum impurities, i.e. local scatterers with 
internal degrees of freedom, into nanofrabricated devices has led to new manifestations of 
Kondo physics, e.g. in quantum dot systems or atoms deposited on metallic surfaces^. In 
these systems many parameters of the impurity, such as its internal spectrum, its coupling 
to the metallic environment and the properties of the latter can be tuned within the ex- 
periment. Investigating the effects of these parameters on observable quantities leads to 
various new questions: originally, most of the theoretical work on the quantum impurity 
problems has neglected the effect of electronic correlations in the host system^. If the lat- 
ter is one-dimensional, however, any interaction leads to non-Fermi liquid behaviour. The 
low energy regime is then described by a Tomonaga Luttinger liquid (TLL) characterized by 
continuously varying exponents of its ground state correlation functions 3 . As a consequence, 
the local density of states vanishes as a power law in such a system. Therefore it is to be 
expected that the critical properties of the impurity will be strongly affected. This problem 
of a quantum impurity coupled to a TLL has been subject of intense studies in recent years - 
both analytically based on field theoretical methods and numerically. These studies indicate 
that depending on the coupling parameters different non-trivial fixed points at intermediate 



or strong coup 
(see e.g. Refs. 



ingcan be realized which determine the observables accessible to experiments 



ffl>. 

This picture of the low energy behaviour of these systems has to be supported by methods 
which do not rely on the analysis at weak coupling, e.g. exact solutions as for the Kondo 
and Anderson impurity problem in a Fermi liquid^. For results which cover the full range of 
experimentally available parameters specific realizations of quantum impurity models have to 
be studied. An approach which allows to make contact to the universal low energy behaviour 
identified using the methods mentioned above is based on integrable lattice models. Starting, 
e.g., from the Bethe ansatz solvable supersymmetric t-J model for electrons on a one- 
dimensional lattice or variations thereof one can consider different representations of the 
graded Lie algebra gl(2\l) for the spectrum of electronic states on bulk and impurity sites 
without destroying integrabilit y 10 ! 11 ' 12 ^ 3 . The physical properties of these impurities can be 
tuned by variation of the representation (which may depend on a continuous parameter for 
gl(2\ 1)) and by a shift in the spectral parameter which directly enters the coupling between 
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the impurity and the host system. 

The resulting Hamiltonian of models constructed along these lines takes a particularly 
simple form when such impurities are combined with open boundarie o 1411 ^ 1 ^ . In this paper 
we shall consider the system introduced in Ref. 3, based on the supersymmetric t—J chain 
with open boundary conditions. The Hamiltonian of this model is 



n 



v j=2 cr 



L 

3=2 



SjSj + i 



HS Z -nN + H b 



(1) 



where V projects out double occupancies on the bulk sites (j = 2 to L + 1) and Sj = 
T^aB c j,a^a/3 c j,/3i n j = c j,a c j,° are the electronic spin and number operators on site j. 
The magnetization of the system is controlled by the field H and we consider the system 
at fixed hole concentration 5=1 — ^2jnj/L (alternatively one may use a grand canonical 
approach to control S by variation of the chemical potential fx). The impurity is added at the 
boundary (site 1) of this chain, its internal spectrum and coupling to the bulk is determined 
by 



At 2 + (a + 2) 2 



(2) 



where Q 2)0 = |0)(a| 2 and Qi,o- = \/a + l|0)(cr|i — 2a\fa |cr)(2| 1 are generalized electron 
annihilation operators for site 1 and 2. The parameter a > labels the four- dimensional 
representation of gl(2\l) used in the construction of (j2j) and controls the internal spectrum 
of the impurity. The terms in (j2J describe exchange and Coulomb interaction between the 
electrons on site 1 and those in the chain as well as a term allowing for the hopping of 
electrons between the bulk and the impurity. Comparing this model with that of the single- 
impurity Anderson model, Vo = 4f2+ ^ +2 - ) 2 can be identified with a hybridization coupling. 
Since we consider an open chain, the parameter t can be either real positive or purely 
imaginary, thereby allowing to cover the entire range — 00 < Vo < 00 for this coupling 
between the bulk and the impurity site. Note that additional parameters can be introduced 
into the model by adding static boundary fields^. For generic parameters, however, the 
1/L contributions of the impurity, boundary and boundary fields to the thermodynamic 
properties are additive in the integrable model. Since the t—J model with open boundaries 



and boundary fields has already been studied in great detail we restrict our analysis to the 
impurity contributions to the magnetization and the susceptibilities. These are functions 
of the bulk density of holes 5, the magnetic field H and the parameters controlling the 
impurity, i.e. a and Vq. Furthermore we consider only the ground state properties of the 
system to avoid the subtleties arising in the analysis of systems with open boundaries at 
finite temperatures^. 

In the following section we identify the configuration of the impurity in the ground state 
of the many particle system as a function of these parameters. Then we study the magnetic 
susceptibility and occupation of the impurity at zero magnetic field where we derive analytic 
expressions near half filling (5 = 0) which are complemented by numerical results for general 
5. Based on these results we identify the relevant energy scales in the system which allow 
for a quantitative study of the properties of the system in an external magnetic field. 



II. FORMATION OF BOUND STATES 



A. Bethe Ansatz Equations 



The spectrum of the model has already been obtained by means of the algebraic Bethe 
ansatz technique 17 . An eigenstate with N e = + electrons is characterized by the roots 
of the Bethe ansatz equations (BAE): 

M s M c 

(e 1 (X j )) 2L = J] e 2 {X j - X k )e 2 (X 3 + \ k ) J] e -i( A i - ^) e -i( A i j = 1; • • • M a> 

Ms [6) 

e a (tf 7 + t)e a (^ -t) = J\ e_i(0 7 - A fe )e_i(tf 7 + A fc ) , 7 = 1, . . . M c . 

fc=i 

Here e n (x) = The M s = L + 1 — JV-j spin rapidities Xj parameterize magnetic 

excitations starting from the completely filled state with maximum polarization (1 electron 
per bulk site, 2 electrons on the impurity site with = L + 1, iV| = 1) and the M c = 
L + 2 — N e charge rapidities $ 7 describe holes added to this state. The phases e a in the 
second set of these equations are due to the presence of the impurity (J2J), similarly as in 
Ref. The energy of the corresponding Bethe state is then given by the expression 

E = V (a + 2) + 2(L - 1) - ^ + - y ) M c + HM S - + 2) - jL. (4) 



The ground state and the low-energy magnetic and charged excitations (spinons and 
holons) of the t—J model with open boundary conditions without boundary fields or an 
impurity - similar to the model with periodic boundary conditions - are described by positive 
rapidities {Ay,i? 7 } solving the BAE (jHJ). Sufficiently strong boundary magnetic fields or 
potentials lead to the formation of boundary bound states in the spectrum of the system. 
In terms of the many-particle Bethe states this is reflected by the appearance of isolated, 
purely imaginary roots to the BA E 20 ' 21 ' 22 . Note that a necessary condition for such roots 
is the existence of singularities in the boundary phase shifts to compensate divergencies 
appearing in the scattering phases for imaginary rapidities. A similar mechanism for the 
creation of bound states has been found to exist in integrable impurity problems^, although 
it has not been studied systematically so far. For the present model with the impurity 
described by Eq. (J2J) bound states appear for a sufficiently strong hybridization Vo- This 
regime is reached by choosing t = ir purely imaginary (r > without loss of generality) . In 
the presence of a particular solution with an imaginary root Im(i? 7o ) > 0, the r.h.s. of the 
second set of equations in (j^J) vanishes in the thermodynamic limit (M s — ► oo with M SjC /L 
kept fixed). Therefore, $ 7 has to be exponentially close to a zero of e a ($ 7 + £)e a (i? 7 — t) = 
e Q +2r('^7)e Q -2T('^7), i-e. $ 70 = — «(a/2±r). This bound state solution appears for r > a/2 = 
t . Increasing r further eventually leads to the appearance of additional imaginary roots 
in the ground state configuration. The analysis of this sequence of bound states suggests a 
division of the V§-a parameter space of the impurity into four regions labeled by an index 
(R) counting the number of bound state solutions present in the ground state of the system 
(see Fig. QJ. R runs from (the region described by the BAE (J3|) ) to III (a region with 
three bound states). 

(I): r <r < (a+l)/2 = n 

As seen above, a complex root $m c = i(r — a/2) appears in the set of the charge 
rapidities {$ 7 } for r > a/2. Explicitly taking into account this root, describing a 
bound state in the holon sector, and rearranging the equations we end up with a set 
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of modified BAE: 

M. 



e i L (\) = e -(2r+i-a)(Aj)e 2r _i- a (A j ) J^e 2 (Aj - X k )e 2 (Xj + X k ) 

k+i 

Mc-l 

x II e -^ X i ~ ^) e -i(^ + #p) , j = l,..., M„ (5) 

13=1 

M s 

e a +2r(^7)e a -2r(^ 7 ) = JJe_i(0 7 - A fc )e_i(^ 7 + X k ) , j = 1, . . . , M c - 1 . 

k=l 

(II): T 1 <T<a/2 + l = r 2 

Increasing r even further, we will generate a new bound state, this time, in the spin 
sector. We associate to this bound state, the complex spin rapidity Xm s = i(r — (a + 
l)/2). Once again, the BAE have to be changed in consequence: 

M a -1 

e i L (^j) =e-(2r-a-3)(Aj)e 2T -i-a(A j ) e 2 (Xj - X k )e 2 (X j + X k ) 



k^j 

M c -1 

X 

/9=1 

M.-l 



JJ e_ 1 (A i -^)e_ 1 (A i + ^), j = 1, . . . , M s - 1, (6) 

/3=1 

Af,-1 

e a+2r (^ 7 ) =e 2r _ 2 _ Q (i9 7 ) j [ e_i(^ 7 - Afc)e_i(i? 7 + A fc ) , 7 = 1, . . . , M c - 1. 

k=l 

From (J3J) we see that the contribution of this 'spinon bound state' to the energy of 
this Bethe state will be 

Eu. = [(t - a/2)(l - t + a/2)]- 1 . (7) 

(III): r > r 2 

For t > a/2 + 1, a third bound state is created. Another purely imaginary root, 
i?Af c -i = *( r — a/2 — 1), will coexist with a set a real rapidities in the charge sector 
and with the root of region (I). The new BAE in this region are: 

Ms-l M c -2 

e? L (A J )= e 2 (A, -A fc )e 2 (A, + A fc ) e_ 1 (A i - ^^(A,- + fy) , 

fc^i ,3=1 

j = l,...,M s -l, (8) 

M s -1 

e a+ 2r(tij) =e 2r _ 2 _ Q (^ 7 ) e_i($ 7 - A fe )e_i($ 7 + A fc ) , 7 = 1, . . . , M c - 2. 

fe=i 



The bound states in this sector, consisting of two charge and one spin rapidity can be 
interpreted as a singlet with vanishing charge bound to the impurity. 

Increasing r beyond a/2 + 1 does not lead to additional bound states as expected for an 
impurity with a single orbital allowing for occupation of at most two charges. Therefore, we 
conclude that the maximum number of bound states allowed to develop themselves upon 
variation of the hybridization at the boundary impurity site is three - two holons and one 
spinon. 



B. Continuum limit - Equations for the densities 

The analysis of the BAE is simplified by doubling of the real roots of the BAE with pos- 
itive and negative ones identified, i.e. X-j = —Xj and $_ 7 = — $ 7 . In the thermodynamic 
limit, the real roots {Xj} ({"$ 7 }) form continuous distributions which are conveniently de- 
scribed in terms of their densities p s {X) (p c (?9) ) . Those densities obey the following coupled, 
but linear, integral equations: 






1 pf\X)+pf\X)\ + -f A A ai L> 

Here we have introduced a y (x) = ^ y 2/l +x 2 and f* k f*g denotes the convolution J_ h dyf(x — 
y)g(y). The boundaries of integration A and B for the spin and charge sector, respectively, 
are fixed by the conditions 



n L (10) 

-B L 

where 9{x) is the Heaviside step function and r& are the thresholds for the appearance 
of bound states identified above. Note that the boundaries of integration A and B are 
completely fixed through (fTUjl by bulk quantities (i.e. total densities of spin a electrons). 
Alternatively, they can be specified in a grand canonical approach by the conjugate poten- 
tials, i.e. the chemical potential p and the magnetic field H. In this case the ground state 
is obtained by filling all modes with negative dressed energies e a (A) and e c (i9) solving the 



integral equations 

e.(A)\ = f -2ir ai (\) + H\ + ( - f^a 2 f^aA # / e s (X) 

with e s (±A) = and e c (±S) = 0. The unpolarized ground state for vanishing magnetic 
field, for example, corresponds to A = oo. 

Since © is a linear system of integral equations, the generic solution for both spin and 
charge densities is of the form 



P = Poo + -j\Pim V + Pb)- 



(12) 



The first term in (jl2j) is the bulk density obtained by solving Q with L = oo. The remaining 
two terms of order 1/L are the contributions due to the presence of the impurity and due 
to the openness of the chain, respectively. 

Depending on the region fixed by the parameters a and t (or r) the driving terms p s , c 
in the integral equations © are obtained from the discrete BAE (jHJ), (JEJ, © or (jHJ) as 



P ( s R) W 



R = 

(A) + a 1+a _ 2r (X) R = I 

— «3+a-2r(A) — a 2T -a-l(A) R = II 

R = III 



(13) 



for the spin-sector and 



(14) 



a a ($ + t) +a a ('&-t) R = (t = irER) 
a a+ 2r($) + a a _2r(tf) R = (r = —it G R) 

flQ+2r(^) + a Q _ 2r (^) i2 = / 
k «a+2r(^) + a2+a-2r(^) # = ^ 

for the charge-sector. The contributions due to the boundaries can be calculated with the 
driving terms pi (A) = 02(A) and ph\'0) = — ai(i?). As a consequence of the decomposition 
()12j) . the bulk, impurity and boundary contributions to any thermodynamic quantity can 
be studied separately. For instance, the ground state energy per site is formally 



+ 



-imp 



+ e b 



L 



(15) 



From this, magnetization, density of electrons and the corresponding susceptibilities are 
obtained by taking the appropriate derivatives with respect to the conjugate fields. The 
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focus of this paper is on the characterization of the impurity. Therefore, we will concentrate 
on the contributions to these quantities which are determined by (j!3|) and (|14|). The finite 
size contributions of the boundaries are of purely geometric nature and have been studied 
by EBler 18 . 



C. The impurity's ground state 

So far, we have distinguished four regions of the a-r parameter space of the impurity, 
characterized by different possible bound state configurations. Each of theses regions is 
described by a different set of BAE. Before studying how the properties of the system are 
affected by the presence of the impurity, the configuration corresponding to the true ground 
state has to be identified. For this, the impurity's contribution to the energy has to be 
computed using the applicable sets of BAE in the vicinity of the transition lines, r = 
The formal expression for this contribution in terms of the densities is 



e 



imp 



A 



(a + 2)V - 2 - 7T / dA pf >(A)ai(A)+ 



-A 



G //^^) + ^-to) + £(t-t 2 )) (16) 
+ # Q / dA (A) + 9(r - n)) - 2//. 

Here, we want to concentrate on the case H = 0: for small magnetic fields (i.e. A ^> 1), a 
convenient treatment of the integral equations © is obtained by rewriting the one for p s as 
follows 

„<*, = gL + ( r A + n * , P <«. + fj^ t (17 , 

1 + a 2 V-oo / 1 + 02 1 + a 2 

Now, at i7 = (A — >• oo), the densities are given in terms of the solution to a scalar integral 

equation for the density pc of charge rapidities: 

B ( 18 ) 



Pc —Pc,H=0 + / ^1 * Pc 
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The driving terms at zero field are 



R = 

G2r-aW + Ga-2rW R = I 

— G2+ a -2rW ~ G2 T ~a-2\^) R = II 

R = III 



(19) 



and 



Pc 



R = 

,t" + G 2r _ a+1 (^) + G a _ 2T+1 (tf)) i2 = J 

G3+«-2t(0) - Gzr-a-!^) R = II 
\, P i ' R = III 

for all four regions. The functions p c have been defined above in ()14|) and 



5 {II) - 

Pc 



(20) 



Gp(x) 



da; . W-g|aV2|) 

— expf— tujx) —-i — -— . 

2tt FV ' 2cosh(u;/2) 



(21) 



Now, we can express the impurity's contribution (|16|) to the energy in terms of the solution 
to the integral equation fT%|) . 

(0) : t = it real and < r < r 

In the absence of bound state, the 1/L correction to the ground state energy reads 

= & + I fj® [p-2TtG l {$)\pf\$) 1 (22) 
where E b = V (a + 2) - 2(p + 1). 

(1) : t < t < n 

As we have seen previously, in region (I) a first bound state is created in the charge 
sector. The energy contribution due to the impurity now becomes: 

-B 



e \mp — Eb — Vt[G2t-q+i(0) + G a -2r+l 



-imp 



\J tofr-^^pPifl+li. (23) 



Here p c has to be evaluated with the appropriate driving term ({2Dj) for region (I). The 
additional chemical potential is due to the charge in the bound state and the terms 
containing the G-functions are the consequence of the rearrangement of the rapidities 
in the spin sector. 
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(II): n < r < r 2 

In region (II), two bound states are possible. Using the appropriate driving terms in 
(|18|) . the 1/L correction to the energy becomes: 

1 /" B 

42 = + 7r[G 3+a -2r(0) + G 2r - a -i(0)] + - / dtf [// - 27rG 1 (^)]p^ ) (^) +H + E Ma - 

* J-B 

(24) 

The extra term, E Mg , is the energy contribution (J7J) of the spin bound state. 

(Ill): t>t 2 

In region (III), €j mp takes the form 

^ = E b + \ j\#[»-27rG l (#))pi ni \#)+2 f , + E Ms . (25) 

The proper ground state configuration is the one which minimizes the impurity contri- 
bution to the energy as given in (|2"2~|) - (j23|) . In Fig. we present numerical data for €i mp , 
for fixed a = 1 and bulk hole concentration 5 = 0.2, as a function of the hybridization Vo 
(note that e^ p can be continued to regions with larger R, describing a configuration where 
an allowed bound state is not occupied). From these numerical data we conclude that - as 
far as the ground state of the impurity is concerned - three intervals in the hybridization 
have to be distinguished (see also Fig. QJ): 

Vo < l/(a + 1) (real t and r = — it G [0, r ] U [r 2 , oo]): for attractive (Vq < 0) and weakly 
repulsive (Vo > 0) hybridization the bound states identified above are not occupied, 
hence the ground state is described by Eqs. (jUJ) with R = 0. 

l/(a + 1) < Vo < A /(2a + 3) (i.e. r < r < rj): in this interval the holon bound state is 
occupied in the ground state configuration is described by with the driving terms 
for R = I. 

Vq > 4/ (2a + 3) (i.e. t\ < r < r 2 ): for strong repulsive coupling between the impurity and 
the host the ground state is obtained using the configuration from R = II, i.e. with 
the holon and the spinon bound state occupied. 

Note that the charge Q = singlet bound state of region III is never present in the ground 
state of the system. 

In the following we continue to label our results by the index R = to III, the relation 
to the physical ground state as a function of Vo, however, is given by the classification above. 
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III. ZERO-FIELD SUSCEPTIBILITY AND OCCUPATION OF THE IMPURITY 



A. Analytical results close to half-filling 



In this section, we will explicitly calculate the magnetization of the impurity in a small 
field. With our parameterization of the BAE's roots, the impurity contribution to the 
magnetization is given by 

M imp = \j\$p¥\$) + \{9{T-T Q ) + 9{T-T 2 ))- 1 - f\\p^\X)-9(r-r 1 ). (26) 

Proceeding as for Eq. (fTTj) . we can rewrite this expression for the impurity's magnetization 
as an integral over the spin density only: 



1 f°° 1 
M imv = - / dA pf >(A) - 9(t - n) + -(9(r - to) + 9(t - r 2 )). 



(27) 



Introducing g(z) = p{ R \A+z) in the integral equation (JUJ), we obtain the following equation 
for the unknown function g: 



POO 

9(z) = p ( s % (A + z) + dz' G 1 (z - z')g(z') 

Jo 

POO pB 

+ / dz' G 1 (2A + z + z')g(z') + / dz' G (z - z' + A)pf\z'). 

JO J-B 



(28) 



For small magnetic fields corresponding to large values of A this equation can be solved by 
iteration using Wiener- Hopf methods. Following Ref. Q, we expand g = g\ + g 2 + where 

POO 

g x {z) = gf\A + z) + C^G (A + z) + / dz' G 1 (z - z') gi (z') (29) 

Jo 

POO POO 

9n(z) = dz' G 1 (2A + z + z')g n „ 1 (z')+ / dz' G x {z - z')g n {z') , n>l. (30) 
Jo Jo 



Here g^ = p^# =0 and we have used that, for A B, 

B pB 

dz' G (A + z- z')pf\z') ^G (A + z) / dz' e™' pf\z') 

B J-B 



(31) 



to define the number C^} = dz' e nz ' p[ R \z') which is given in terms of pi alone. 

The leading behaviour of the impurity magnetization for large A is now obtained from 
()29|) : using the results (|A5|) and ()A6|) from Appendix El we find 



■^J-imp 



V2 



ire 



C { a} T + 2 COS 



■ tt(2t— q) 1 



Cg ) +2cos(^); 



R = 
} +1/2 R = I 
-1/2 R = II 



(32) 
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Here the dependence of M imp on the density of electrons in the host is completely given 
through the constants C ajT . The zero field limit of these quantities is given in terms of the 
solution to the integral equation (fT%|) for pi . In general this equation has to be solved 
numerically. For B <C 1, |2r — a\ (i.e. close to half-filling and the impurity sufficiently far 
away from the threshold for the holon bound state), however, it can be solved by iteration. 
Doing so, we obtain at first order in B the following expressions for C ar : 



Ci R } = 2B I 



■la 



R = 

l+a-2-rl*-' » R = I 

G 3+ a-2r(0) - G^ a ^(0) R = II 



7r(a 2 — At 2 ) 

■k(J-At 2 ) + ^2r+l-a(0) + Gl+a-2r(0) 

4(l+q) 
7r(2+a-2T)(Q+2r) 



(33) 



At the threshold, r = r = a/2, the leading contribution to Ci°| , Ca,l at small hole 
concentration are ±1 independent of 5. 

Finally, the boundaries of integration, A and B, have to be expressed in terms of phys- 
ical quantities, namely the concentration 5 = M c /L of holes (doping) in the bulk and the 
magnetic field using the relations ()1U|) and Again, restricting ourselves to the regime 

close to half filling we have tt5 = IB In 2. To express A in terms of the magnetic field one 
has to enforce e s (A) = 0. A Wiener- Hopf analysis of the integral equations give. 



g 18 



nA = - HH/H ) + 



4 In if' 



Hn 



xf(2^/ej(2n - C) 



(34) 



with C = f_ B d$ e^Sc^d), here e c is the dressed energy of the holes given by (fTTj). Close to 
half-filling C < 2n and H ps (2n) 3/2 / \fe. Using (jHSJ) and (JH1J) the low-field magnetization 
is found to be linear in H . Thus we obtain for the impurity contribution to the zero field 
magnetic susceptibility close to half filling close to half filling 5 1, |2r — a|: 



(0): 



(I): 



a5 



Ximp 



2^ 



Ximp 



7l5 



vr 2 (a 2 -4r 2 )ln2 



+ 0(S 2 



(35) 



4a 



cos 



2 In 2 \n(a 2 -At 2 ) 
7t(2t — a) 



+ G 2t+ l-a (0) + G 1+Q _ 2T (0) 



(36) 



+ 0(S 2 ) 
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(II): 



Ximp 



7l5 



27T 2 
+ COS 



4(1 + a) 



2 In 2 V vr(2 + a - 2r) (a + 2r) 
T(2r-a)', + ( ^ 



G*3+q-2t(0) — G*2r-a-l(0) 



(37) 



Finally, we can compute the number of electrons on the impurity which is obtained from 
(ITU1) to be 



D 



imp 



2 - 



B 



(38) 



B 



Again this expression can be calculated at zero field, close to half-filling and away from the 
holon bound state threshold (B <C 1, |2r — a|) giving 

' 2- S, .,J a ,.o, R=0 
L-^rl — +^ (V.,._, _„(0) G',_„_-v(0)] # = / 



D 



imp 



2 In 2 V 7r(a 2 -4r 2 ) 



+ G 2r+1 _ a (0) + G' 1+Q _ 2T (0)) 



(39) 



1 - 



4(l+a) 



2 In 2 \n(2+a-2T)(a+2T) 



G 3+ a-2r(0)-G aT _ a _i(0)) fl = // 



B. Numerical results for arbitrary doping 

For finite density of holes (B finite) the integral equations (fT%|) and (fTTjl for the charge 
components of the densities and dressed energies have to be solved numerically In Figures 01 
HI and El we present results of this analysis for the zero field susceptibility and occupation of 
the impurity for different doping as a function of the hybridization Vq. Note that for general 
filling factors, the constant C in the definition (J3l|) of H is no longer negligible, thus 



Ximp(H = 0) 



1 

2n(2n - C) 



a 



(o) 



7r(2r— a) > 
2 - 



CS + 2 cos( 



cS ) + 2cos(^^; 



(40) 



R = 
R = I 
R = II 

First notice that the a = case where the doubly occupied state decouples from the 
remaining three electronic impurity states |0), \a) is special: here, for Vo < 1, the impurity 
is doubly occupied D imp = 2. Therefore the impurity does not contribute to the magnetic 
susceptibility of the system, i.e. Ximp = 0. For Vo > 1 (regions (I) and (II)), the occupation is 
less than 2 due to the filled holon bound state and there is a small finite impurity contribution 
to the susceptibility. 

For a > 0, the total number of electrons present at the impurity can become larger than 2 
as a consequence of the attractive interaction between the impurity and the bulk electrons for 



14 



V < 0. Note, however, that the difference D imp \ Vo =-oc, — D imp \y 0=+oo approaches 2, showing 
the depletion of the impurity orbital at large positive Vq. For small hole concentration 5 we 
had found above that the impurity occupation changes at the threshold for the formation 
of the holon bound state Vq ^ l/(a + 1). For finite 5 this resonance moves to smaller values 
of Vo, determined by the condition that the (real) impurity parameter t is of the order of 
the boundary B which is given by the bulk hole concentration through Eq. Q. Below this 
value of the hybridization, i.e. for t > B the impurity is doubly occupied and essentially 
decoupled from the host. 

The same shift is observed in the resonance of the impurity contribution to the zero-field 
susceptibility which moves from the threshold for the holon bound state towards smaller 
values of Vq as the hole concentration is increased. This resonance is the response of the un- 
paired electron on the impurity site which appears when D imp « 1. Although still limited by 
fluctuations in the impurity's occupation the susceptibility at the resonance grows strongly 
with 5. This is shown in Figure El where the maximum of Ximp{H = 0) as a function of Vq is 
given for different values of the impurity parameter a as a function of the hole concentration 
in the bulk t-J chain. For 5 — *> the susceptibility approaches that of the bulk system while 
it diverges for 5 — > 1, i.e. vanishing bulk density of electrons. In this limit the remaining 
electron on the impurity site is essentially an uncoupled local moment. 

IV. IMPURITY IN A MAGNETIC FIELD 

As we have seen above the magnetic behaviour of the impurity is most interesting in the 
regime of weak coupling Vq parameterized by real values of t. This parameter introduces 
a scale where the response of the impurity to external fields is expected to change. In the 
present problem where the impurity site is coupled to the t-J chain with separated spin and 
charge degrees of freedom this response will depend on the relation of this scale to the typical 
energies in these sectors as well. First, using the relation ([34)1 . we obtain the Kondo field 
determining the scale where the susceptibility of an impurity coupled to the spin degrees of 
freedom changes as 



for large t corresponding to Vo <C 1. Note that the dependence of Hk on the hybridization 
is exponential and not a power-law as found for a Kondo impurity in a TLL4»^. This reflects 




(41) 
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the absence of backscattering in the integrable impurity model considered here. 

On the other hand, as seen from the BAE (|5]L the primary effect of the impurity (J2J) is to 
induce a phase shift in the charge sector. Hence, the hybridization has to exceed a minimal 
value below which the impurity is doubly occupied and will not contribute to the magnetic 
response of the system (this is different from the Anderson impurity model considered by 
Bortz et al. where the local orbital can be populated by at most one electron^). Finally, for 
the impurity to be visible in the magnetic response, spin and charge sectors have to be cou- 
pled at the scale introduced by the hybridization. This coupling is determined by the relative 
size of the impurity parameter t and that of the host parameters A, B. Note that the latter 
are functions of the magnetic field: we have already used that A — > oo while B approaches 
a finite constant for H — > 0. On the other hand, the system becomes completely polarized 
for H = H sat = 4cos 2 (7r5/2). In this limit A -> A sat = \^/A/H sat - 1 = §tan(7T<?/2) 
while B — > oo. Therefore, depending on the relative size of A and B two scenarios can be 
distinguished for fixed hole concentration 5 in the bulk: 

(i) For sufficiently weak hybridization, corresponding to large t such that t = B > A 
at some value H = Hb of the magnetic field, the impurity is decoupled from the (charge 
degrees of the) host at the Kondo scale (JHJ). In this case the impurity contribution Ximp to 
the susceptibility will exhibit a resonance at H = H B > H K while being suppressed below 
H B . 

(ii) If the impurity and the host are already coupled at the Kondo scale, i.e. t = A < B, 
Kondo like behaviour of the susceptibility can be observed where the susceptibility ap- 
proaches some non-universal finite value as H — > 0. 

Finally, for Vq large enough for Hr > H sat corresponding to t < A sat (or t purely 
imaginary) no resonance is present in the susceptibility due to the finite bandwidth of the 
lattice model. 

In Figure the regions corresponding to these scenarios are shown for various values of 
the hole concentration 5. By numerical integration of the Bethe ansatz equations (jUJ) the 
occupation of the impurity in these regions is determined. Fluctuations in the densities 
are measured by the impurity contribution to the magnetic susceptibility Ximp and charge 
compressibility K imp = dD imp /dfi of the system. Just as the corresponding bulk quantities 
these thermodynamic coefficients are conveniently computed based on their representation 
in terms of the dressed charge matrix (see Appendix [Bj) • In Fig. |H1 our numerical results 
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on these quantities as a function of the hybridization and the magnetic field are shown for 
hole concentration 5 = 0.8. For intermediate values of the hybridization 0.1 ^ Vo ^ 0.2 
the coupling of the impurity to the holon excitations is effective. For small magnetic fields 
below the Kondo scale Hk the impurity contribution to the susceptibility takes a non- 
universal value xo characteristic for the strong coupling regime of an Anderson impurity. 
Above Hk the field dependence is that of a local moment with logarithmic deviations from 
full polarization (see Figure EJ). The emergence of universal Kondo like behaviour Ximp — 
f(H/H K )/2nH K for smaller hybridization V < 0.1 is suppressed by the decoupling of the 
impurity from the host. Here the double occupancy of the impurity for small fields and 
the formation of a local moment with an impurity occupation close to 1 appears for fields 
H > Hb 3> Hk are clearly visible. In the vicinity of the transition between these regions the 
occupation and magnetization of the impurity fluctuate strongly, as shown by the resonances 
in the susceptibilities (Fig. |HJ). Finally, for larger values of V > 0.2, the occupation of the 
impurity decreases well below 1 and the susceptibilities are approximately constant over the 
entire range of the external field < H < H sat . 

V. SUMMARY & CONCLUSION 

We have studied the properties of an Anderson-type impurity embedded into a super- 
symmetric t-J chain with open boundaries. Within this model the hybridization coupling 
Vo of the impurity to the holons in the ID host can be tuned freely while preserving the 
integrability. Upon variation of Vo the nature of the many-particle ground state changes 
due to the appearance of a sequence of bound states in the holon and spinon sectors. From 
our analysis of the contribution of the impurity to the electronic density and magnetization 
of the system and the corresponding susceptibilities at vanishing magnetic field we have 
identified two regimes: for attractive or weakly repulsive hybridization the impurity is dou- 
bly occupied with vanishing fluctuations. Here the impurity and the host are effectively 
decoupled due to a mismatch of Vq and the relevant scales in the holon sector. For strong 
hybridization one holon and one spinon bound state are occupied. Between these regions 
the susceptibility exhibits a sharp resonance which diverges as the density of bulk electrons 
tends to 0. In a magnetic field exceeding the characteristic scale for excitations in the spinon 
sector, H > Hk, a third region appears for intermediate values of the hybridization: this 



17 



regime is characterized by the formation of a local moment with small corrections to full 
polarization. In this range of Vq the field dependence of the susceptibility approaches the 
characteristic Kondo scaling behaviour until the decoupling of the impurity from the holon 
sector sets in. 

Note that the zero field analysis of Section IIHI can be extended to small magnetic fields. 
This will lead to additional logarithmic corrections to the thermodynamic quantities. In the 
response of the Hamiltonian they will be hidden by the ones due to the presence of open 
boundaries 18 . 

The model considered in this paper can be generalized in a number of ways. As discussed 
in the introduction, the boundary impurity can be combined with a local potential or mag- 
netic field which by itself will generate a sequence of bound states in the holon and spinon 
sector. Fine-tuning the parameters describing the impurity and those of the boundary field, 
part of the impurity spectrum can be projected out from the Hilbert space 25 . This proce- 
dure leads, e.g., to exactly solvable models for a true Kondo spin S or an impurity coupled 
only to one spin channel of the host electrons (see also Refs. fl^jlfiL While the projection 
onto a subset of impurity states is well understood on the level of the Hamiltonian and its 
construction it is an open problem, how the corresponding spectra (which will differ for 
different choices of the projected subset) are to be computed. In particular the role of the 
bound states in these sectors needs further investigation which is left for future work. 
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APPENDIX A: SOLUTION OF THE INTEGRAL EQUATIONS (|5U|> 

The integral equations (J30|) are of Wiener- Hopf (WH) type 




(Al) 
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and can be solved using standard techniques based on the factorization of the Fourier trans- 
formed kernel Gi(u) 

[l-diuj)]- 1 = G + (u)G-(uj) , UmG ± {uj) = l (A2) 

into functions G^iu) which are analytic for Im(uj) > (< 0), respectively. For the t—J 
model these techniques have been applied befor e - 1 ^ and the factorization of the kernel is 
known to be 



G-(u) = G+{-u) = ^ 2 \ (Z:) ■ (A3) 



In Eqs. f!29|) and (|30j) for n = 1 and 2 three different driving terms go need to be considered: 
(a) The case go(z) = G (A + z) already appears in the calculation of bulk quantities such 
as the dressed energies for the t—J model. Following Refs. we find for g 

9 +{ u )=iG + {u;)G-{-in)-^-. (A4) 

uj + in 

Using the explicit expressions (|A3jl we obtain 

2 -kA 




dz g a (z) = g+(u = 0) = \ —e~* A , (A5) 



en 



which will be necessary to compute the impurity's magnetization. 

(b) The second type of driving term, according to ([Uj) , is go(z) = Gp(z+A)+G-.p(z+A). 
The analysis of the WH equation is completely analogous to the first case and we find 

g+(u) = 2iG + (uj)G-(-iir)-^cosA = 2 cos A 9 +(uj). (A6) 

(c) Finally, the driving term in the equation (j30|) for the sub-leadin g co ntribution, g2, 
is proportional to go(z) = J °° dz' G±(2A + z + z')g a {z'). Following Ref. bja, we perform a 
Laplace tansform of go{z) to obtain: 

1 f 00 

g (z) ~ — J dx e- 2Ax e~W x gt(ix) (x + . . .) , (A7) 

where we have used the asymptotic expansion of the function G\(z) ~ 1/47T2; 2 + 0(z~ 4 ). 
Now the solution of the Wiener-Hopf equation is given by 

~ G + (uj)— dx e- 2Ax (x + . . .) G+ ( ta; )g°M (A8) 
47r Jq u; + ix 
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The presence of the rapidly decaying factor exp(— 2 Ax) (remember that A ^> 1 at small 
field) in the integrand suggests the following expansion around x = 



-wA 



g+(ix)G + (ix) (x + . . .) ~ 2^=x + 0(x 2 ). (A9) 

\/ en 



(A10) 



From this expression we obtain 

APPENDIX B: EXPRESSIONS FOR THE SUSCEPTIBILITIES 

For the numerical analysis of the susceptibilities in a finite magnetic field it is convenient 
to use their expressions in terms of the so called dressed charge matrix 27 ' 28,29 . This quantity 
is as the solution of the Bethe ansatz integral equations 



UA) UA) \ ( 10\ / -J* A a 2 J% ai \ I £ SS (A) 




* 



(Bl) 



(B2) 



J V 01 / I fi a i / 

The magnetic susceptibility at fixed p is obtained by first taking the derivatives of the 
magnetization as obtained from (fTUj) with respect to the boundaries of integration A and B. 
Then, starting from the conditions e s (A) = = e c (B), one obtains two equations for their 
derivatives OA/ OH and OB /OH which are then solved in terms of £ at the boundaries of 
integration 

7 Z ss Z cs \ ( £ SS (A) £ CS (A) 

Z sc Zee 

This approach leads to the following expression for the bulk magnetic susceptibility at fixed 
chemical potential 2 ^ 

Xtulkl = ^ f ^- 2Z ^ 2 + ( Z -- 2Z -) 2 \ (B3 ) 

Alt y V c V s J 

where and v s = e' s (A)/7ip s (A), v c = e' c (B) / n p c (B) are the spinon and holon Fermi ve- 
locities, respectively. The impurity contribution to the magnetic susceptibility is obtained 
analogously starting from (}2T)j) as 10 

, = 1 f {Z cc -2Z sc ff c + {Z cs -2Z ss y fs \ 
M 47T V v c v s J 

with f s = p s R \A)/ Ps (A) and f c = p c R) (B) / p c (B) . 
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In this paper we work at fixed (bulk) hole concentration 5 = d$ Poo,cC#)- Using 
85/ OH together with the expression of the chemical potential entering Eqs. in terms of 
the independent variables H and 5 the bulk susceptibility is found to be (see also Refs. 21 2^) 



_ (det Zf 1 
Xbulk ^-^^v c Zl + v s Zl- (B5) 

Again, it is straightforward to compute the impurity contribution to the susceptibility within 

this approach giving 

Ximpls = tPl — ; { fsZ cc {Z ss — -Z cs ) — f c Z cs {Z sc — -Zee) \ ■ (B6) 

n V c /j cs + V s^cc I 1 1 ) 

To measure the valence fluctuations on the impurity site one has to consider the charge 
compressibility. Again, the bulk and impurity contributions to k = dN e /dfi are conveniently 
expressed in terms of the dressed charge matri x 1 ^ 2 ^ 2 ?! 2 ? 



ry 2 ry 2 
K'bulk \h r 



TTVc TTV S , R7 s 

7 2 y 2 ( B7 ) 

I ^CC p ^CS r 

^imp ij Jc i J s ■ 

IX V r 7TVa 
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FIG. 1: Spectrum of impurity bound states in the Vq-ol parameter space. 
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FIG. 2: Impurity contribution to the ground state energy for different bound state configurations 
as a function of the hybridization Vq for bulk hole concentration of 5 = 0.2. The impurity parameter 
is fixed to a = 1. Note that the Bethe ansatz for the configuration R = I and R = II gives identical 
results for Ei mp at negative Vq. 
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FIG. 3: Zero field susceptibility Ximp (upper panel) and electron number Di mp on the impurity 
site (lower panel) for bulk doping 5 = 0.2 as a function of Vq. 



26 




FIG. 4: Same as Fig. for 8 = 0.5. 
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FIG. 5: Same as Fig. for 8 = 0.8. 
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FIG. 6: Maximum value of Ximp{H = 0) as a function of 8. For small 6 this maximum appears 
close to the threshold for the creation of the holon bound state Vq = l/(a + 1). 
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FIG. 7: Relevant energy scales for the impurity problem as a function of the hybridization Vq for 
a = 1: for fixed hole condensation 5 (e.g. full lines for 5 = 0.8) the impurity is decoupled from the 
host for magnetic fields below the left line (shaded area). Above this threshold (or at sufficiently 
large Vb) the Kondo scale Hk (right branch of full line) becomes visible in the impurity's response 
to the external field. The dotted lines indicate the continuation of the Kondo scale Hk (|41jl into 
the decoupled region. 



30 




0.1 0.2 0.3 



FIG. 8: (a) Susceptibility Ximp-, (b) electronic occupation Di mp , and (c) charge compressibility 
K-imp (normalized to its bulk value) of the impurity as a function of the hybridization and magnetic 
field for 5 = 0.8 and a = 1 (phase boundaries from Fig. are superimposed). Note the logarithmic 
scale used for the shading of the susceptibilities. 
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FIG. 9: Magnetic field dependence of the impurity susceptibility for 5 = 0.8, a = 1 and various 
values of the hybridization: for Vq ^ 0.1 the transition between the strong coupling behaviour 
for small fields and the formation of a local moment above H ~ Hk is clearly seen. At smaller 
values of Vq the susceptibility is strongly suppressed due to the decoupling of the doubly occupied 
impurity from the host (data for Vq = 0.05 are enhanced by a factor of 10). 
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